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ABSTRACT 

We report new estimates of the time delays in the quadruple gravitationally lensed quasar 
PGl 115+080, obtained from the monitoring data in filter R with the 1.5-m telescope at the 
Maidanak Mountain (Uzbekistan, Central Asia) in 2004-2006. The time delays are 16.4 days 
between images C and B, and 12 days between C and A1+A2, with image C being leading 
for both pairs. The only known estimates of the time delays in PGl 115 are those based on 
observations by Schechter et al. (1997) - 23.7 and 9.4 days between images C and B, C and 
A1+A2, respectively, as calculated by Schechter et al., and 25 and 13.3 days as revised by 
Barkana (1997) for the same image components with the use of another method. The new 
values of time delays in PG 1115+080 may be expected to provide larger estimates of the 
Hubble constant thus decreasing a diversity between the Hq estimates taken from gravitation- 
ally lensed quasars and with other methods. 
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1 INTRODUCTION 

As was first suggested by Refsdal (1964), gravitationally lensed 
quasars can potentially provide an estimate of the Hubble constant 
Hq independent of any intermediate distance ladder. This can be 
made from measurements of the time delays between the quasar 
intrinsic brightness variations seen in different quasar images. The 
value of Hq can be obtained then (within the adopted cosmolog- 
ical model) from the observed geometry of the system, with the 
known lens and source redshifts, and with the use of physically 
validated model of mass distribution in the lensing galaxy. Since 
a phenomenon of gravitational lensing is controlled by the surface 
density of the total matter (dark plus luminous), it provides a unique 
possibility both to determine the value of Hq and to probe the dark 
matter abundance in lensing galaxies and along the light paths in 
the medium between the quasar and observer. 

By now the time delays have been measured in about 20 grav- 
itationally lensed quasars resulting in the values of Hg, different 
for different quasars, while remaining noticeably less than the most 
recent estimate of Hg obtained in the HST Hubble Constant Key 
Project with the use of Cepheids - Hq = 72 ± 8 km s"' Mpc"' 
(Freedman et al. 2001). This discrepancy is large enough and, if 
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the Hubble Constant is really a universal constant, needs to be ex- 
plained. A detailed analysis of the problem of divergent Hq esti- 
mates inherent in the time delay method and the ways to solve it can 
be found, e.g., in Keeton & Kochanek (1997); Kochanek (2002); 
Kochanek & Schechter (2004); Schechter (2005). 

The quadruply lensed quasars, and PGl 1 15+080 in particular, 
are known to better suit for determining the Hq value as compared 
to the two-image lenses since they provide more observational con- 
straints to fit the lens model. The PGl 1 15 source quasar with a red- 
shift of zs = 1.722 is lensed by a galaxy with za = 0.31 (Henry 
& Heasley 1986; Christian et al. 1987; Tonry 1998), which forms 
four quasar images, with an image pair Al and A2 bracketing the 
critical curve very close to each other. It is the second gravitation- 
ally lensed quasar discovered over a quarter of century ago, at first 
as a tripple quasar (Weymann et al. 1980). Hege et al. (1980) were 
the first to resolve the brightest image component into two images 
separated by 0.48 arcsec. Further observations (Young et al. 1981; 
Vanderriest et al. 1986; Christian et al. 1987; Kristian et al. 1993) 
have provided positions of quasar images and information about 
the lensing object, which allowed to build a model of the system 
(e.g., Keeton, Kochanek & Seljak 1997). In particular, Keeton & 
Kochanek (1997) have shown that the observed quasar image posi- 
tions and fluxes and the galaxy position can be fit well by an ellip- 
soidal galaxy with an external shear rather than by a just ellipsoidal 
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galaxy or a circular galaxy with an external shear. They noted that 
a group of nearby galaxies detected by Young et al. (1981) could 
provide the needed external shear. 

The problem of determining the Hubble constant from the 
time delay lenses is known to suffer from the so-called central con- 
centration degeneracy, which means that, given the measured time 
delay values, the estimates of the Hubble constant turn out to be 
strongly model-dependent. In particular, models with more cen- 
trally concentrated mass distribution (lower dark matter content) 
provide higher values of /fo> more consistent with the results of the 
local Hq measurements than those with lower mass concentration 
towards the centre (more dark matter). 

The time delays in PG 1 1 15-1-080 were determined for the first 
time by Schechter et al. (1997) to be 23.7 ± 3.4 days between B 
and C, and 9.4 ± 3.4 days between A1-I-A2 and C (image C is 
leading). Barkana (1997) re-analyzed their data using another al- 
gorithm and reported 25^3 g days for the time delay between B and 
C. There are also estimates of the time delays between images Al 
and A2 made from the Chandra and XMM-Newton X-ray Observa- 
tories data (Dai X, et al. 2001, Chartas et al. 2004). As is predicted 
by all lens models, it does not exceed a small fraction of the day, 
and equals 0.16 ± 0.02 days (Chandra data) or 0.149 ± 0.006 days 
(XMM-Newton data). Determination of the time delays has gener- 
ated a flow of models for the system, (Schechter et al. 1997; Keeton 
& Kochanek 1997; Courbin et al. 1997; Impey et al. 1998; Saha & 
Williams 2001; Kochanek, Keeton & McLeod 2001; Zhao & Pronk 
2001; Chiba 2002; Treu & Koopmans 2002; Yoo et al. 2005, 2006; 
Pooley et al. 2006; Miranda & Jetzer 2007), all illustrating how 
strongly the estimated value of Hq depends on the adopted mass 
profiles of the lens galaxy for the given values of time delays. 



2 OBSERVATIONS 

Observations of PG 1115+080 have been carried out at the 1.5- 
meter telescope of the high-altitude Maidanak Observatory (Cen- 
tral Asia, Uzbekistan) with a scientific BroCam CCD camera. A 
SITe ST005A 2030 x 800 chip provided an image scale of 0.26"/pix 
at the f/8 focal plane. The CCD images were usually taken in series 
consisting of 2 to 10 frames for the R filter and of 2 to 6 frames for 
V and /. To provide higher photometric accuracy, we averaged the 
values of magnitudes estimated from individual frames. The seeing 
varied from 0.75" to 1.3" (the FWHM of images of the reference 
stars B and C according to designation by Vanderriest et al. 1986). 
The analysis of photometry shows no significant dependence of the 
photometry errors on seeing, excepting the FWHM noticeably ex- 
ceeding 1.3". Occasional frames with such values of the FWHM 
were excluded from processing. 

The algorithm of photometric image processing is similar 
to that applied to photometry of Q2237+0305 and described in 
great details in (Vakulik et al. 2004). The light curves of PG 
1115+080A1,A2,B,C in filter/? for three observational seasons in 
2004-2006 are shown in Fig. 1. The data sets consist of 23 nights 
in 2004, 27 in 2005, and 24 in 2006. The data demonstrate notice- 
able variations of the quasar luminosity, with the total amplitude 
reaching almost 0.4 mag during a year (between April 2004 and 
May 2005). The time delays between the light curves of the C and 
B, C and Al (or A2) images can be easily seen from a simple vi- 
sual inspection of the R light curves. Thus, our photometry in filter 
R has immediately demonstrated applicability of the obtained light 
curves to determine the time delays. As compared to the data used 
by Schechter et al. (1997) and Barkana (1997), we were lucky to 



detect the quasar brightness variations with an amplitude of almost 
a factor of three larger, and with rather well-sampled data points 
within every season of observations. In addition, the accuracy of 
our photometry has made it possible to confidently detect flux vari- 
ations with an amplitude as small as 0.06 mag, that can be seen in 
the data of 2004. 

Also, recently published by Morgan et al. (2008) observations 
of PGl 115-1-080 in filter R during the same time periods in 2004- 
2006 should be mentioned here. We have made use of their pho- 
tometry presented in their table 3 to compare with our light curves. 
Variations of the quasar brightness which allowed us to determine 
the time delays are seen in their A1-I-A2 light curve quite well, but 
become undetectable in the B and C light curves because of a much 
larger scatter of the data points. 



3 CALCULATION OF THE TIME DELAYS 

The ideology of methods to determine the time delays between two 
image components from their light curves is simple enough and ob- 
vious. A common feature of all known methods of time delay mea- 
surements is the use, in one way or another, of the cross-correlation 
maximum or mutual dispersion minimum criteria, while the meth- 
ods may differ in the algorithms of the initial data interpolation. 
Analysis of the light curves of quasar images in pairs can be also 
applied when a lens consists of more than two images. But however, 
another approach seems to be more promising in this case. The 
model source light curve can be determined from a joint analysis 
of light curves of all image components. The individual time delays 
of the components are then determined with respect to this model 
source light curve jointly from a corresponding system of equa- 
tions. In some cases, this approach allows systematic variations in 
the light curves to be revealed and taken into account, such as those 
caused by, e.g., microlensing. A similar approach described earlier 
by Press et al. (1992) and Rybicki & Press (1992) and used later 
by Barkana (1997) was also applied to determine the time delays in 
Q2237-H0305 (Vakulik et al. 2006). 

As was noted above, a necessity to properly interpolate the 
unevenly sampled data points in the light curves under considera- 
tion is one of the main technical problems in determining the time 
delays. A variety of interpolating functions and algorithms is used, 
such as polynomials of various power, Legendre polynomials, etc. 
In some cases the low-frequency splines provide good results. But 
unfortunately, all these interpolation procedures do not contain any 
physical meaning. Meanwhile, the Fourier spectra of quasar vari- 
ability are known to be rapidly decreasing functions, (de Vries et 
al. 2005), and this is naturally explained by the finite physical sizes 
of quasars. In particular, the quasar size is known to play a role 
of smoothing factor in microlensing light curves, as is clearly seen 
from simulations. Therefore, we tried to find such an algorithm to 
represent the source quasar light curve, which would take into ac- 
count the expected frequency characteristics of the quasar variabil- 
ity and allow a relevant smoothing of the observational data. To 
do this, it was natural to address an expression that is well known 
in optics, radio engineering and theory of information as the sam- 
pling theorem. According to this theorem, a signal with a bounded 
spectrum can be represented accurately enough by a function: 

f{t)= ^ a(kAt)si?icin(t - kAt)/At], (I) 

Here, the function sinc(x) is specified as sinc(x) = sin(x)/x, and 
a(kAt) are samples of f(t) taken at a mesh with a step At, which is 
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Figure 1. The light curves of PG 1 115+080 Al, A2, B, C from observations in filter R with the 1.5-m telescope of the Maidanak Observatory in 2004, 2005 
and 2006. 



determined by a boundary frequency Qj,,,/ of the function f{t) spec- 
trum: At = l/(2Q.t„d)- Since, as was noted above, Fourier spectra of 
quasar variability are rapidly decreasing functions, we may apply 
the samphng theorem to reproduce the quasar light curve. We used 
this approach earlier to analyze statistics of microlensing brightness 
variations in Q2237+0305 (Minakov et al. 2008). 

The same algorithm was used in this work to represent the 
model for the source light curve. Since the time delay between the 
light curves of images Al and A2 is very short (Dai et al. 2001, 
Chartas et al. 2004), their fluxes were summed to form a single 
curve, which we will call the A light curve. Thus, we may write the 
following functional for three light curves: 



1 V> V' K(f.) + dmj - f(t„ At, Tj)f 



(2) 



where mjitj) are the data points in the light curve of the j-th image 
at the time moments f,, dnij and Tj are the shifts of correspond- 
ing light curves in stellar magnitude and in time. A' is a number 
of points in the light curves, (Tj(f,) are the photometry errors, and 
finally, /(f,, Af, Tj) is an approximating function (1). 

We adopted dnio = and Tq = in our calculations, that is, 
we fitted the light curves of the two other images to the A light 
curve, and thus, dmi and dm2 are the magnitude differences A-B 
and A-C, respectively. At given values of Tj and T2, we minimize 
<I>(Af, Ti, T2) in rfm, and in coefficients a(kAt) of the sampling func- 
tion. The values of minimum of 0(Af, ti , T2) were being looked for 
at a rectangular mesh Ti,T2 with a step of 0.5 days in preliminary 
calculations, and of 0.2 days at a final stage. The values of Ti , T2 
corresponding to the minimal value of <I)(Af, ti, T2), were adopted 



as the estimates of the time delays Tba and tac- The time delay 
Tfic is not an independent quantity in our method, and can be deter- 
mined as a linear combination Tgc = t^a + Tac- 

Having calculated the time delays, we analyzed deviations 6 j 
of light curves of each image from the approximating function: 



6j(ti) = mj(ti) + dmj - fit,. At, Tj). 



(3) 



We revealed a linear trend in deviations of the A light curve, 
and small parabolic trends for images B and C. We interpreted 
these trends as the effects of microlensing, which are expected to 
be rather slow and weak. We then subtracted a half of these trends 
from the initial light curves of the components, obtained new esti- 
mates of the time delays, and analyzed the residual trends again. 
This procedure continued iteratively until the trends became in- 
significant. The final estimates of Tj and T2 are obtained with the 
linear trend 0.0105 mag/year subtracted from the image A light 
curve, and with the parabolic trend with the amplitude of +0.01 
mag subtracted from the B and C light curves. Thus, our esti- 
mates of the time delays are Tba = 4.4;^2^, Tac = i'^-Ol'^o' 
Tbc = 16.4+24 days, with the relationship TacItba, more consistent 
with that determined by Barkana (1997) than by Schechter et al. 
(1997). 

To evaluate the errors of estimating the time delays and re- 
liability of our estimates, we fulfilled a numerical simulation. We 
selected a function /(/,, Ar, Tj) used to approximate our light curves 
with At = 0.08 yr, as a model source light curve. The model light 
curves of the components were obtained by shifting /(r,, At,Tj) by 
the proper time delays Ti,T2 and magnitude differences, and by 
adding some random quantities to imitate the photometry errors. 
The estimates of these errors were obtained from the analysis of 
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Figure 2. The data points of the A, B and C hght curves superposed with each other after shifting by the corresponding magnitude differences and time delays 

tbc = 16.4+2J obtained in this work. The paraineter At of the approximating function (shown in solid curve) is 0.12 years. 



deviations of the observed data points from the approximating func- 
tion resulting in the values of standard deviations 0.008, 0.016 and 
0.011 mag for images A, B and C, respectively. Since the method 
we used might be susceptible to the mutual locations of data points 
in the actual light curves, we selected the model samples exactly 
at the same time moments as in the actual light curves. We sim- 
ulated two cases: t^a = 14.3 and t^c = 9.4 days as determined 
by Schechter et al. (1997), and Tba = 4.4 and Tac = 12 days (our 
result). We simulated 2000 random light curves synthesized as de- 
scribed above, and calculated the resulting time delays using the 
procedure, which was exactly the same as in the analysis of the 
actual light curves. Admitting that we could be mistaken in selec- 
tion of the Nyquist interval, we fulfilled the model calculations with 
a more low-frequency function of the source brightness variations 
(Af = 0.12 years). No systematic biases larger than 0.3 days in 
the estimates of simulated time delays Ti and T2 were revealed in 
both cases. The results of simulations were used to estimate the 
95-percent confidence intervals. 

It is interesting to note that using only the data of 2004, where 
a small-amplitude turn-over in the light curves is detected, we ob- 
tained Tba = 5.0, Tac = 9.4, and Tbc = 14.4 days, consistent with 
the estimates obtained from the whole data set. But however, sim- 
ulation of errors for only this time interval demonstrates notice- 
ably larger uncertainties, as compared to those calculated from the 
whole light curve. 

In Fig. 2, the light curves of images A, B and C shifted by 
the corresponding time delays and reduced to image A in magni- 
tude are shown for the approximating function parameter Af = 0.12 
years. As is seen, the data points for all the three images are very 
well consistent with each other and with the approximating curve, 
(a behavior of the approximating function within the gaps between 
three seasons of observations should be ignored). 

So we obtained the time delay values, differing noticeably 
from those by Schechter et al. (1997) and Barkana (1997), which 
are used in a variety of models of many authors to derive the Hub- 
ble constant value. Calculation of new lens model or recalculation 
of some most popular ones to derive the new estimate of the Hub- 



ble constant would be well beyond the scope of the present short 
communication. We would just note that the new values of time 
delays reported in this work must result in higher values of Hq as 
compared to those obtained with the previous time delay values for 
PGll 15-1-080, thus decreasing a well-known diversity between the 
time-delay method of determining the Hubble constant and other 
methods. 
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